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Context. The solar irradiance is known to change on time scales of minutes to decades, and it is suspected that its substantial fluctua- 
tions are partially responsible for climate variations. 

Aims. We are developing a solar atmosphere code that allows the physical modeling of the entire solar spectrum composed of quiet 
Sun and active regions. This code is a tool for modeling the variability of the solar irradiance and understanding its influence on Earth. 
Methods. We exploit further development of the radiative transfer code COSI that now incorporates the calculation of molecular lines. 
We validated COSI under the conditions of local thermodynamic equilibrium (LTE) against the synthetic spectra calculated with the 
ATLAS code. The synthetic solar spectra were also calculated in non-local thermodynamic equilibrium (NLTE) and compared to the 
available measured spectra. In doing so we have defined the main problems of the modeling, e.g., the lack of opacity in the UV part 
of the spectrum and the inconsistency in the calculations of the visible continuum level, and we describe a solution to these problems. 
Results. The improved version of COSI allows us to reach good agreement between the calculated and observed solar spectra as 
measured by SOLSTICE and SIM onboard the SORCE satellite and ATLAS 3 mission operated from the Space Shuttle. We find that 
NLTE effects are very important for the modeling of the solar spectrum even in the visual part of the spectrum and for its variability 
over the entire solar spectrum. In addition to the strong effect on the UV part of the spectrum, NLTE effects influence the concentration 
of the negative ion of hydrogen, which results in a significant change of the visible continuum level and the irradiance variability. 

Key words. Line: formation - atomic data - Molecular processes - Sun: atmosphere - Sun: UV radiation - Radiative transfer 



1. Introduction 

The solar radiation is the main source of the input of energy 
to the terrestrial atmosphere, so that it determines Earth's ther- 
mal balance and climate. Although it has been known since 
1978 that the solar irradiance is not constant b ut instead varies 
on scales from several m inutes to decades (cf. lFrohlichll2005l 
■ iKri vova & Sol anki 2008), the influence of this variability on the 
climate is not yet fully understood. Nowadays several datasets 
for the past spectral solar irradian ce (SSI) based on the different 
recon struction approaches (e.g. iLean et ail 120051; iKrivova et alj 
2009) and satellite measurements are available. However, the re- 
maining disagreements between these data lead to different at- 
mospheric response s when they are used in the climate models 
(Sha piro et alJl2010h . The task of constructing a self-consistent 
physical model in order to reconstruct the past solar spectral ir- 
radiance (SSI) remains of high importance. 

Modern reconstructions of the SSI are based on the assump- 
tion that the irradiance changes are det ermined by the evolu- 
tion of the solar su r face magnetic field (Foukal & Leanl fl988; 
IKrivova et al.1 120031; iDomingo et ail |2009). Areas of the solar 
disk are associated to several components (e.g. quiet Sun, bright 
network, plage, and sunspot) according to the measured surface 
magnetic field and the contrast of the features. These compo- 
nents are represented by correspondi ng atmosphere structures 
(cf. lKuruczll99Tl;lFontenla et al.|[T9 99). The SSI is calculated by 
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weighting the irradiance from each model with the correspond- 
ing filling factor. 

The calculation of the emergent solar radiation, even from 
the atmosphere with known thermal structure, is a very so- 
phisticated problem because consistent models have to account 
for the NLTE effects in the solar atmosphere. As the impor- 
tance of these effects has become clear over the past sev- 
eral decades, several numerical codes have been developed. 
One of the first NLTE-codes, LINEAR, was published by 
lAuer et a! .] (119721). who used a complete li nearization method 
developed by lAuer & Mihaiasl (fl969i [l 970). Later, the MULTI 
code was published bv ICarlssonl (Il986l) . This code i s base d on 
the linearization technique deve loped by IScharmeil (Il98ll) and 
IScharmer & Carlssonl d!985albl) . More recently, the RH code, 
which is based on the M ALI (Multi-level Approximate La mbda 
Iteration) formalis m of iRvbicki & Hummed dl99ll 1 19921) has 
been developed bv lUitenbroekl d200 ll) . 

Combining th e ra diative transfer code by 

lHamann & Schmutzl d !987l) ; ISchmutz et alj dl989l) and spectral 
synthe sis code SYNSPEC bv lHubenvl dl981l) . lHaberreiter et all 
(2008) has developed the ID spherical symmetrical COde 
for Solar Irradiance (COSI). The NLTE opacity distribution 
function (ODF) concept, which was implemented in this code, 
allows an indirect account for the NLTE effects in several mil- 
lion lines. This makes COSI especially suitable for calculating 
the overall energy distribution in the solar spectrum. In this 
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paper we introduce a new version of COSI (version 2), describe 
the main modifications, and present first results. 

In Sect. 12 we describe the datasets of the measured spectral 
irradiance, which were used for comparison with our calculated 
spectrum. In Sect. [3] we introduce the integration of molecular 
lines into COSI and show that it can solve the discrepancies to 
observations and to other codes. In Sect. |4]we present the re- 
sulting solar spectrum from NLTE calculations. Due to missing 
opacity in the UV, the flux appears to be significantly h igher than 
measured dBusa et alj200Ulshort & Hauschildtl2 009). We solve 
this problem by introducing additional opacity to the ODF for se- 
lected spectral ranges (Sect. 14. il l, while the problem of the NLTE 
visible continuum due to the deviations in the concentration of 
the hydrogen negative ion is addressed in Sect. 14.21 In Sect. [5] we 
present the synthetic spectra of the active regions and its impli- 
cations for the solar variability study. Finally, we summarize the 
main results in Sect. [6] 



2. Measured solar spectra irradiance 

We compared our calculated solar spectrum with several 
available datasets. We used the observations taken with 
SOLSTICE (SOLar-STe llar Irradiance Comparison Experiment, 
iMcClintock et al.l d2005l)) up to 320 nm and the Solar Irradiance 
Monitor (SIM) dHarder et al.ll2005l) fr om 320 nm onw ard instru- 
ments onboard the SORCE satellite (Rottman 2005) obtained 
during the 23rd solar cycle minimum (hereafter SORCE mea- 
surements). For the comparisons shown in this paper, we used 
the average of the observations from 21 April 2008 to 28 April 
2008. 

Our calculated spectrum was converted with a 1 nm boxcar 
profile for comparison with the SOLSTICE measurements and 
trapezoidal profile for comparison with the SIM measurements. 
The SIM resolution strongly depends on wavelength (FWHM is 
about 1.5 nm for the 310 nm and about 22 nm for the 800 nm) 
so the parameters of this trapezoi dal pro f ile are also wavelength 
dependent and were provided bv lHarden d2009l) . 

We also u sed the SOLar SPECtral Irradiance Measurements 
(SOLSPEC) dThuillier et al.f2 004) during the ATLAS 3 mission 
in November 1994. For this comparison we convolved the cal- 
culated spectrum with a Gaussian with FWHM = 0.6 nm. 
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3. Molecules in the COSI code 

COSI simultaneously solves the equations of statistical equilib- 
rium and radiative transfer in ID spherical geometry. All spec- 
tral lines in COSI can be divided into two groups. The first group 
comprises about one thousand lines that are the most prominent 
in the solar spectrum. These lines are explicitly treated in NLTE. 
The seco nd group contain s several million background lines pro- 
vided by iKuruczl ([2006) and calculated under the assumption 
that their upper and lower levels are populated in LTE relative to 
the explicit NLTE levels. These background lines are taken into 
account in the spectral synthesis part of the code but also affect 
NLTE calculations via the ODF. The ODF is iteratively recal- 
culated until it becomes self-consistent with pop ulations of the 
NLTE levels dHaberreiterl2006tlHaberreiter et al.l2008h . This al- 
lows us to indirectly account for the NLTE effects in the back- 
ground lines. 

Molec ular lines were not incl uded in the previous version 
of COSI ( Hab erreiter et al.l 120081) . However, molecular lines 
play an important role in the formation of the solar spectrum. In 
some spectral regions they are even the dominant opacity source. 
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Fig. 1. The ratio R of the number concentration of carbon and 
oxygen that are not attached to the molecules to their total ele- 
ment amount (upper panels) and the CN and CO concentrations 
as a function of height (lower panels) for three atmospheric mod- 
els: FALA99 (dashed curves), FALC99 (continuous curves), and 
FALP99 (dotted curves). The zero-height depth point in all mod- 
els is defined as the radius at which the continuum optical depth 
at 5000A is equal to one. 



Therefore, a considerable amount of opacity was missing. This 
led to a discrepancy with observations. In this version of COSI 
molecular lines were taken into account. This requires comput- 
ing molecular concentrations, preparing the molecular line lists, 
and incorporating molecular opacities and emissivities into the 
code. 



3. 1 . Chemical equilibrium 

Under the assumption of the instantaneous chemical equilib- 
rium, the concentration of molecule AB is connected with the 
con centrations o f atoms A and B by the Guldberg-Waage law 
(cf. lTaturifl966h : 



«A»B 
«AB 



( 2nmkT \ V2 Q A Q B 
\ h 2 J Qab 



Kab(T) 



kTexp(-D /kT), (1) 



where K^s is the temperature dependent equilibrium constant, Q 
internal partition function for the atoms A and B, and molecule 
AB, m is the reduced molecular mass, and Do the dissociation 
energy. The contribution to the opacity from the transitions be- 
tween lower level i and upper level j of molecule AB is given 
by 
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where g, is statistical weight of the lower level, B's are the 
Einstein coefficients, and <p(v) is the line profile normalized to 
one. 

There are several available temperature polynomial approx- 
imations for the atomic and molecular part ition functions, as 
well as for th e equi librium constants Kar (c f .lTatumlll966l:lTsuiil 
[19731: llrwrn1[T98l ISauval & Tatuml fl98l iRossi et al] 119851 
Although molecular partition functions can contain significant 
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errors due to the inaccurate values of the ene rgy levels and even 
due to unknown electronic states (llrwinl 19811) . this does not sig- 
nificantly affect the accuracy of molecular opacity calculations. 
This insensitivity exists because the overall molecular concen- 
tration is proportional to the partition function (see Eq.Q]), while 
the opacity in a given line is inversely proportional to it (see 
Eq. |2j. Therefore it is crucial to either calculate equilibrium con- 
stants directly from partition functions or use the same source for 
bot h approximat ions . 

Tsuii d 19731) shows that abundances of hydrogen, carbon, 
nitrogen, and oxygen can be determined quite accurately with- 
out taking into account their dependencies on other elements. 
Because we are only interested in molecular bands that can sig- 
nificantly contribute to the opacity over a broad spectral range, 
the system of equations like Eq. (Q]) were solved only for these 
four main elements and the diatomic molecules built from them. 

Because a significant fraction of atoms can be associated 
with molecules the chemical equilibrium calculation also af- 
fects the atomic lines. In Fig. Q] we present the changes in 
the carbon and oxygen concentrations due to the association 
with molecules, together with CN and CO concent rations for 
three atmospheric models by iFontenla et alj ([1999): the rela- 
tively cold supergranular cell center model FALA (hereafter 
FALA99), the averaged quiet Sun model FALC (hereafter 
FALC99), and the relatively warm bright network model FALP 
(hereafter FALP99). Both molecular concentrations and devia- 
tions in atomic concentration show strong temperature sensitivi- 
ties, and this is important for assessing of the solar spectral vari- 
ability. 

3.2. Main molecular bands in the solar spectrum 

The line list for the spect rum synthesis was compiled using 
the molecular databases of iKuruca (1 19931) and Solar Radiatio n 
Physical Modeling (S RPM) database (e . g.lFontenla et al.ll2009l) . 
which is based on iGrav & Corballvl (11994 and fflTRAN. 
Wavelengths and line strengths of the most significant lines of 
the CN violet system and CH G band w e re calculated based o n 
the molecular con s tants b v lKxuppl(u974l) . lKnowles et all (11988). 
and lWallace et al.l d!999l) . T he OH and CH conti nuous opacities 
were calculated according to lKurucz et al.l (Q987) under the LTE 
assumption. 

In Fig. |2] we present the part of the solar spectrum where 
the LTE calculations with the previous version of COSI, which 
did not account for molecular lines, showed significant devi- 
ations from the calculations with th e LTE radiative transfer 
code ATLAS 12 carried out bv lKuruczl(l200l . Both calculations 
used t he same atmosphere structure and abundances by iKuruczl 
(119911) (hereafter K91). One can see that introducing the molec- 
ular lines into COSI quite strongly affects the spectrum and sig- 
nificantly diminishes disagreements with the ATLAS 12 code. 
The most prominent features are the CH G band around 430 nm, 
CN violet G band around 380 nm, CN, NH, and OH bands be- 
tween 300 and 350 nm. The remaining differences can possibly 
be attributed to the different photoionization cross sections im- 
plemented in COSI and ATLAS12. 

In Fig. [3] we compare the spectra calculated by both codes 
with the SORCE measurements. One can see that both codes 
predict the correct value of the continuum in the red part of 
the visible spectrum (600 nm) but fail to reproduce the molec- 
ular bands and UV part of the spectrum. The LTE spectra 
calculate d with the FALC 99 and two sets of abundances by 
iGrevesse & Anders! d 199 lh (hereafter G9 1 ) and bv lAsplund et al.1 
(120051) (hereafter A05) are presented in Fig. [4] The calculations 
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Fig. 2. LTE calculation of molecular lines with COSI. Upper 
panel: Synthetic solar spectrum calculated with (solid line) 
and without (dashed line) molecular lines. Lower panel: 
Corresponding ratios to the irradiance calculated by the 
ATLAS 12 code. All spectra are averaged with a 1-nm boxcar. 



2.5 F 



B 
d 

CS) 

B 

O 

d 

53 

fc! 




200 300 400 500 
Wavelength [nm] 



600 



Fig. 3. The irradiance calculated by COSI in LTE (thin solid line) 
and the ATLAS 12 (dashed line), using K91 atmosphere struc- 
ture and abundances. The spectrum observed by SORCE is given 
by the thick solid line. All spectra are convolved with the instru- 
ment profile of SORCE. 



with the G91 abundances predict the correct value of the visi- 
ble and IR continua. Furthermore, they give quite a good fit of 
molecular bands (see right panel of the picture where the calcu- 
lation of CN violet system and CN G band are shown in more 
detail). However the calculated irradiance in the UV part is again 
much higher than observed for both abundance sets. 

The metal abundances by A05 are significantly lower (up to 
about 2 times) than in G9 1 . The lower metal abundance leads to a 
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significantly lower electron density, consequently to a lower con- 
centration of negative hydrogen ion. The latter is the ma in source 
of the continuum opacity in the solar atmosphere (cf. iMihalasI 
Il978l p. 102). Therefore the application of the A05 abundances 
results in a lower continuum opacity for the visible and IR wave- 
lengths and as a consequence in a higher irradiance. This effect is 
especially prominent at longer wavelengths (from about 400 nm) 
where the number of strong lines is relatively low and the opac- 
ity of negative hydrogen dominates all other sources of opacity. 
In turn, the lower metal abundance given by A05 leads to lower 
line opacities and photoionization opacities in the UV. As a re- 
sult, the A05 abundances also lead to a higher irradiance in the 
UV than the G91 abundances. The FALC99 model has so far 
mainly been used with G91 abundances. 



4. NLTE calculations 

Most of the solar radiation emerges from regions of the atmo- 
sphere with a steep temperature gradient. Therefore a signifi- 
cant part of the escaping photons is not in the thermodynamic 
equilibrium with the surrounding medium. This implies that a 
self-consistent radiative transfer model has to account for NLTE 
effects. It is well-known that these effects are extremely impor- 
tant for the strong lines and the UV radiation. In this section we 
present NLTE calculations of the solar spectrum with COSI and 
show that the NLTE effects are also important for the overall 
energy distribution in the solar spectrum, including the visible 
and infrared wavelength ranges. We want to emphasize here that 
the NLTE calculations presented in this paper were performed 
using temp erature and pressure p rofiles of the solar atmosphere 
obtained bv lFontenla et al.1 ([1999). These profiles in turn depend 
on assumptions, so all effects presented in this paper are rela- 
tive to the assumed model of the solar atmosphere (see also the 
discussion at the end of Sect. [43}. 

In the left panel of Fig. [5] the NLTE calculations of the solar 
spectrum with the atmosphere structure FALC99 and the abun- 
dance sets A05 and G91 are compared to the observed solar 
spectrum. The right panel of Fig. [5]shows the ratios of the NLTE 
to the LTE calculations. The most prominent overall NLTE ef- 
fects are the strong increase in the irradiance in the UV and mild 
decrease of the irradiance in the red part of the spectrum. The 
first effect is discussed in detail in Sect. 14.11 and the second one 
in Sect. 1431 



4.1. Calculating of the UV radiation 
4.1.1. NLTE opacity effects 

Independent of the solar atmosphere model and abundances, the 
LTE calculations overestimate the UV irradiance up to about 300 
nm (see Figs.[3]and[4]i. This problem is even more severe in the 
NLTE calculations as they predict a UV irradiance that is higher 
than in the case of the LTE calculations (see left panel of Fig. [5]). 

The effect of the increase in the f ar UV flux in the NLTE cal - 
culations is actually well known (cf. IShort & Hausch ildt 2009). 
The UV radiation incident from the higher and hotter parts of the 
solar atmosphere onto the photosphere causes a stronger ioniza- 
tion of iron and other metals relative to the LTE case. This de- 
creases the populations of the neutral atoms and consequently 
weakens the strength of the corresponding line and photoioniza- 
tion cross section, which are the main sources of the opacity in 
the far-UV. Therefore it leads to a decrease in the opacity and 
consequently an increase in the flux in the wavelengths up to 
about 250 nm (red threshold of the photoionization from the Mg 



I 2 level). In the longer wavelengths the line and the photoion- 
ization opacities of the ionized metals become more important 
so the increase in their populations results in lower irradiance. 

The effect of the increased UV flux in the NLTE calculations 
can be decreased significantly by the use of the opacity distri- 
bution function (see Fig. [6j, which allows incorporation of the 
opacity from the backgroun d lines into the solu tion of statisti- 
cal equilibrium equations (IHaberreiter et al.ll2008l) . This back- 
ground opacity diminishes the amount of the penetrating UV 
radiation and decreases the degree of metal ionization. It leads 
to the relative decrease in the spectral irradiance up to 250 nm 
and an increase at higher wavelengths (in agreement with the 
above discussion). The complicated wavelength dependence of 
the spectrum alterations caused by the influence of the ODF can 
be explained by the large number of the strong metal lines and 
photoionization thresholds. Overall, introducing the ODF causes 
the redistribution of the solar irradiance over the entire spectrum 
as it decreases the far-UV part and increases the near-UV and 
visible. 

Although the concept of the ODF allows the agreement be- 
tween calculations and observati ons to be imporved, t he ODF 
based on the Kurucz linelist (IHaberreiter et alJ I2008I) cannot 
solve the discrepancy with measurements. The most probable 
source of the problem is the inaccuracy of the atomic and molec- 
ular line list and possible missing continuum opacity. The pro- 
duction of a complete line list is an extremely laborious task. 
Considerable progress has been reached during past decades 
thanks to the effort of the OPAC ITY project dSeatonl I2005I) 
and, in particular, of iKuruczl d!993l) in calculating atomic data. 
Concerning the problem of the UV blanketing, it is important to 
keep in mind that about 99 % (!) of the atomic and molecular lines 
are predicetd theo retically and o nly 1% of all lines are observed 
in the laboratory (lKuruczll2005l) . This leads to potential errors 
in the total opacity as lines could still be missing while wave- 
lengths and oscillator strengths of the predicted lines could be 
inaccurate. An incompleteness of the atomic data is especially 
significant in the UV where the immense number of lines form 
the UV line haze. 

The consistent account for the line blanketing effect is 
very important for th e overall absolute flux distribution espe- 



cially for the UV (c f. jCollet et al 



2005; Avrett & Loeser 2008; 



iFontenla et aT]|2009l) . lBusa et all d200ll) used version 2.2 of the 
MULTI code, which UV opacity were contributed only from 
several hundred calculated in the NLTE lines, bound-free and 
free-free processes. They showed that NLTE calculations with- 
out proper account for the line blanketing effect can lead to over- 
estimating of the UV flux by six orders of magnitude in the case 
of cold metal-rich stars. To compensate for this, they multiplied 
the continuum opacity by a wavelength-dependent factor and 
gave an explicit expression for the fact or parameterization. A 
similar approach is used bvlBrulsl d 1993) to ca lculate the NLTE 
populations of Ni L lShort & H auschildti ( 20091) show that the use 
of different line lists for the line blanketing calculations can lead 
to significantly different spectral irradiance distributions. They 
also show that the problem of the excess of the UV flux in the 
300-420 nm spectral region can be solved by increasing of the 
continuum absorption coefficient. 

4.1 .2. Additional opacities in COSI 

COSI t akes the opacity from t he several million lines into ac- 
count (IHaberreiter et alJ I2008I) . However, the significant dis- 
agreement between the calculated and measured solar UV flux 
as discussed above indicates that the computations still miss an 
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Fig. 4. Spectra calculated by COSI in LTE using the FALC99 atmosphere structure and A05 and G91 abundances, vs. SORCE 
measurements. All spectra are convolved with the instrument profile of SORCE. 
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Fig. 5. Left panel: as the left panel in Fig.|4]but calculated in NLTE. Right panel: The ratio between NLTE and LTE emergent flux. 
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Fig. 6. Left panel: NLTE calculations with and without iterated ODF, compared to the solar irradiance as observed by SORCE. 
Right panel: Ratio between the emergent flux from NLTE calculated with and without iterated ODF. All calculations are done with 
the FALC atmosphere model and G91 abundances. 



essential part of the opacity. To compensate for this missing 
opacity. lHaberreiter et alJ J2008) calculated the UV part of the 
solar spectrum with artificially increased Doppler broadening. 
This approach proves itself successful, but it cannot be used to 
fit high-resol ution spectra. Ther efor e, here we have extended the 
approach by iBusa et alJ (1200 lb and IShort & Hauschildtl {2009) 



and included additional opacity in our NLTE-ODF scheme, as 
described below. 

COSI overestimates the irradiance, i.e., shows an opacity de- 
ficiency in the spectral region between 160 nm and 320 nm. 
At shorter wavelengths, the opacity is dominated by continuum 
photoionization, while at longer wavelengths tabulated atomic 
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and molecular lines are sufficient to reproduce the observa- 
tions. In this wavelength region, we grouped the solar spec- 
trum into 1 nm intervals (corresponding to the resolution of the 
SOLSTICE/SORCE spectrum). For each interval we empirically 
determined the coefficient f c {A) by which the continuum opacity 
has to be multiplied to reproduce the SOLSTICE measurements. 
The additional emission coefficient was calculated under the as- 
sumption that the missing opacity source has thermal structure, 
so its source function is equal to the Planck function. Because 
the correction of the ODF affects the statistical equilibrium of 
the populations, we iteratively solve the statistical equilibrium 
equations and the factors f c {A) using newly updated ODFs un- 
til we reached a self-consistent solution. Thus, the f c {A) factors 
were found iteratively. We stopped the iteration when the devi- 
ations to the observed spectrum was smaller than 10%. In the 
considered region there were very few wavelength points where 
the irradiance was underestimated. We did not apply any correc- 
tion to these points. This underestimation can be explained by 
the inaccuracy of the line list because it can be caused by the 
small wavelength inaccuracy of a few strong lines. 

The wavelength dependence of f c (A) is shown in Fig. [7] for 
the calculations with the FALC99 atmosphere model and G91 
and A05 sets of abundances. One can see that f c (A) is very close 
to 1 (which means little additional opacity) for wavelengths 
longer than 280 nm. There are also several prominent peaks 
that indicate a strong lack of opacity at specific wavelengths. 
Although different abundances lead to a change in the peak am- 
plitudes, the overall profile of wavelength dependence stays the 
same for both sets of abundances, G91 and A05, respectively. 
The averaged line opacity is higher by several orders of magni- 
tude than the continuum opacity, so the ratio of the additional 
opacity, which was needed to achieve agreement with the ob- 
served solar spectrum, to the total opacity that was already in- 
cluded into COSI is for most wavelength below 5 % (see lower 
panel of Fig. [7]). Therefore, the correction of the overestimation 
of the emergent UV flux can be achieved by a moderate modifi- 
cation of the input line list. 

Calculations with the G91 abundances lead to a lower con- 
tinuum than calculations with A05 abundances. Therefore the 
flux overestimation is stronger in the case of calculations with 
A05 abundances. Consequently one has to use larger correction 
factors / c (/l) for A05 abundances. The correction factor strongly 
depends on the chosen line list. We tested ou r calculations usini 
the V ienna Atomic Line Database (VALD) dKupka et alj|199' 
2000) and retrieving the lines with known radiative damping 
parameter. Because thi s database contains fewer lines then pro- 
vided bv lKuruca d2006l) . a larger factor f c (A) is needed to achieve 
agreement with the observations 

4.2. Calculation of the visible and IR radiation 

NLTE effects decrease the visible and near IR irradiance by 
about 10 % (see Fig. [5]), because they influence the concentra- 
tion of electrons and the negative hydrogen ions that determines 
the continuum opacity, as well as the continuum source function. 



4.2.1. NLTE effects in electron and negative hydrogen 
concentrations 

Departures from LTE in the electron and hydrogen negative ion 
concentrations are presented in Fig. [8] for the calculations with 
FALC99 atmosphere structure and A05 abundances. The upper 
panel illustrates the main electron sources in the solar atmo- 
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the G91 and A05 abundances. 



sphere. In the lower part of the photosphere and in the chro- 
mosphere, hydrogen is ionized and provides most of the elec- 
trons; however, in the photosphere, which is responsible for the 
formation of the visible and near-IR irradiance, the temperature 
is too low for ionizing hydrogen. Therefore, here the electrons 
are mainly provided by the metals with low ionization potential, 
such as iron, silicon, and magnesium. 

As one can see from the second upper panel of Fig. [8] the 
ionization of silicon and iron is strongly affected by the NLTE 
effects. This is basically the same effect of NLTE "overioniza- 
tion" as already discussed in Sect. 14.11 where the main emphasis 
was on the deviations in the concentration of the neutral atoms. 
Applying of the ODF significantly decreases the NLTE effects 
on the electron concentration as it strongly increases the UV 
opacity. In the last part of Fig. [8] the contribution functions for 
the continuum radiation are plotted for different /i-positions. The 
contribution fu nctions sho w where the emergent intensity origi- 
nates (see, e.g. lGravlll992l p. 151), so one can see that, although 
the continuum is mainly formed in the lower part of the photo- 
sphere where the NLTE effects in the electron concentration are 
not so prominent, they still can be important, especially for the 
radiation emitted close to the solar limb. 

The photoionization cross section of negative hydrogen has 
its threshold at abo ut 16500 A, a nd it reaches its maximum at 
about 8500 A (cf. lMihalaslll978l p. 102). Although the photo- 
sphere is optically thick for the UV and radiation in the strong 
atomic or molecular lines, the visible and IR photons can easily 
escape, even from the lower part of the photosphere. (The op- 
tical depth of the zero point in the continuum radiation is about 
one, depending on the wavelengths.) This means that all over the 
photosphere there is a lack of photons that are able to photoion- 
ize the negative hydrogen. Thus, in the NLTE case, the radia- 
tive destruction rate of FT is decreased, which in turn increases 
the concentration of the negative hydrogen ion compared to the 
LTE equilibrium. This is illustrated in the third panel of Fig. [8] 
The effect is especially strong in the layers where the contin- 
uum radiation is formed (see the contribution functions in the 
lower panel of Fig. [8). Moreover, in contrast to effects on the 
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lines) continuum radiation and temperature on depth dependence 
(solid curve) for the FALC99 atmosphere model. The contribu- 
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tions are carried out for the FALC99 atmosphere model and A05 
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electron concentration, applying the ODF does not decrease the 
NLTE deviations of the negative hydrogen concentration since 
UV radiation does not make any significant contribution to its 
ionizations. The effe ct of negative hydroge n underionization was 
firstly pointed out bv lVernazza et alJ dl98lh . who found that pho- 



tospheric departure coefficients of the negative hydrogen (ratios 
between NLTE and LTE concentrations) are greater than one. 

Thus, compared to the LTE case, NLTE ionization of hy- 
drogen and metals result in higher electron concentration. The 
NLTE ionization of negative hydrogen also leads to the increase 
in its concentrations. This enlarges the continuum opacity so that 
continuum radiation comes from higher and cooler layers of the 
photosphere. In addition, there is also a decrease in the NLTE 
continuum source function, which contributes to the decrease in 
the continuum emission. 

The effects described above significantly decrease the level 
of the emergent flux in the visible and near IR. Although the 
NLTE treatment is physically more consistent than LTE, it leads 
to a severe deviation from the observed solar spectrum (see 
Fig.0). 

4.2.2. Collisions with negative hydrogen 

One possible cause of the inability of the NLTE calculations to 
reproduce the measurements is the inadequate treatment of col- 
lisional rates for the negative hydrogen ion. The present version 
of COSI accounts for the electron detachment through collisions 
with electrons and neutral hyd rogen, as well as charge neutral- 
ization with protons as given bv lLambert & Pagell(ll968l) . We do 
not treat collisions of H with heavy elements, which in princi- 
ple could lead to a strong underestimation of the collisional rates. 
We found that increasing the collisional rates involving FL by a 
factor of ten solves the disagreement between the measured and 
calculated continuum flux for the models using A05 abundances. 
However, when using increased collisional rates, the molecular 
bands (especially CH G band) appear to be weaker than the ob- 
served ones. To make the molecular lines consistent with ob- 
servations we changed the A05 carbon and nitrogen abundances 
and used them as given by G91. This allows us to reach a very 
good agreement with measurements (see Sect. 14.31 ). 

Let us emphasize that an increase in the collisional rates for 
H cannot help adjust the NLTE calculations performed with 
G91 abundances, because the considered increase in collisions 
can only enforce an LTE population ratio between negative hy- 
drogen and H 1 1 level. The latter is, however, in an NLTE regime 
mainly thanks to the NLTE "overionization". As the calcula- 
tions with G91 abundances give good agreement in a purely LTE 
regime (see Fig. @), the considered increase in the collisional 
rates for H~ is not enough to enforce full LTE agreement. An 
NLTE treatment of the solar atmosphere yields population num- 
bers of negative hydrogen that strongly deviate from LTE. 



4.2.3. Change in abundances 

The effect of using of different abundances on the emergent syn- 
thetic spectrum is illustrated in Fig. [9] One can see that the abun- 
dance change of elements like iron, magnesium, and silicium 
strongly affects the UV radiation (mainly due to the photoion- 
ization opacities), but also the visible and infrared continuum, 
because these elements are the main donors of the photospheric 
electrons. In contrast, a change in nitrogen and oxygen abun- 
dances has no effect on the continuum and only affects the spec- 
trum via the molecular bands. The carbon abundance slightly 
influences the continuum, but the main effect of carbon in the 
visible spectral range comes from molecular systems. The CN 
violet system at about 380 nm and CH G band at about 430 nm 
are the most prominent features. Thus, the carbon and nitrogen 
abundance changes from the A05 to G91 values almost do not 
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affect the continuum level (see Fig. |9]l and lead to good agree- 
ment with the measurements (see Fig. ITOb . Therefore, we can 
conclude that A05 abundances should be given preference for 
calculating the NLTE continuum, except for carbon and nitro- 
gen, for which the G91 abundances lead to better reproduction 
of the molecular bands. The favored combination of abundances 
leads to good overall agreement with the observed solar spec- 
trum. 

4.3. Comparison with measurements and discussion 

In Fig. [TO] we present the solar spectrum calculated with the 
FALC99 atmosphere model compared with the SORCE mea- 
surements. The agreement in the 160-320 nm spectral region 
is automatically good (see Sect. 14.1.21 ). At shorter wavelengths, 
the observed flux is determined by several strong emission lines, 
which are treated in LTE in the current version of COSI. This 
leads to a strong underestimation of the irradiance. In a future 
version of COSI, these lines will be treated in NLTE. The Ly a 
line has already been computed in NLTE, and we obtain an over- 
estimation of its flux. In a forthcomin g paper we will investigate 
it in more detail (ISchoell et alj|2010t) . 

The detailed comparison of our synthetic spectrum with 
SOLSPEC measurements is given in Figs. [13] and [14] in the 
Online Material. One can see that starting from 160 nm two 
spectra are in remarkable agreement with each other. Although 
several strong emission lines are currently not properly cal- 
culated in shorter wavelengths, the continuum level there is 
also consistent with SOLSPEC measurements. This can only be 
achieved with NLTE calculations. 

Thus the assumption of enhanced negative hydrogen colli- 
sions and the use of the combined G91 and A05 abundance set 
provide us with the possibility of successfully modeling the en- 
tire solar spectrum. Another possible cause of the problem that 
synthetic NLTE models have difficulty fitting the observed solar 
spectrum is that the FALC99 atmosphere model was developed 
to fit the visible and I R continuum with th e radiative transfer 
code PANDORA (see iFontenla et al] 1 19991) . which uses a dif- 
ferent approximate approach to treat the NLTE effects, differ- 
ent sets of the atomic and molecular data a s well as different 
sets of abundances. The lFontenla e t al. ( 1999) atmosphere struc- 
tures have been use d successfully in vari ous applications (cf. 
iPenza et"aL 2004aB IVitas & Vincel 120051) in which, however, 
the continuum flux was calculated in LTE. The quest to con- 
struct a self-consistent atmosphere model that fully accounts for 
all NLTE effects will be addressed in a future investigation. 

5. Irradiance variations 

It is generally accepted that most of the solar variability is intro- 
duced by the competition between the solar flux decrease due to 
the dark sunspots and a flux increase due to the bright magnetic 
netwo rk and plage (e.g. IWillson & Hudsod[l9?llKrivova eTall 
120031) . It is therefore important to consistently calculate the solar 
irradiance from these active compon ents of the solar atmos phere. 

For our calculations we used the lFontenla et alJ (Il999h atmo- 
sphere model FALC99 for the quiet Sun, model FALF99 for the 
bright network, model FALP99 for plage, and model FALS99 for 
the sunspot. Both NLTE and LTE spectra of the quiet Sun, active 
network, plage, and sunspot are presented in Fig.QT] The bright 
network spectrum is hardly distinguishable from the quiet Sun 
spectra. The NLTE spectra produce a lower visible and UV flux 
and, accordingly, lower total solar irradiance (TSI). This effect 
for the quiet Sun irradiance was already discussed in Sect. 14.21 



Table 1. TSI differences with FALC model per 1 % area (in 
W/m 2 ). 





FALF 


FALP 


FALS 


LTEG91 


0.28 


1.48 


-10.7 


LTE A05 


0.25 


1.27 


-10.82 


NLTE G91 


0.18 


0.76 


-11.02 


NLTE comb. 


0.15 


0.99 


-11.05 



Another interesting detail is the decrease in the contrast be- 
tween the quiet Sun, plage, and bright network in the NLTE cal- 
culations. The main reason for this is that the NLTE effects that 
decrease the continuum source function are stronger in the hotter 
FALP model as it corresponds to the lower particle density (and 
consequently lower collisional rates). 

Figure[T2]shows the flux differences between the active com- 
ponents of the solar atmosphere and the quiet Sun calculated in 
LTE with A05 abundances, and NLTE with enhanced collisions 
and combined A05 and G91 abundances (see Sect.©. The most 
prominent peak at about 3890 A corresponds to the CN violet 
system. Although the differences in the CH G band around 4300 
A are also clearly visible, they are less pronounced than in the 
CN violet system. This is caused mainly by the differences in 
the CH and CN dissociation energies (3.465 eV and 7.76 eV), 
so the CN concentration is more sensitive to temperature varia- 
tions (see Sect. 13.11 ). Furthermore, the contrast between the ac- 
tive components and the quiet Sun is significantly decreased in 
the NLTE calculations. As a consequence, the NLTE calcula- 
tions reduce the solar variability in the visible and IR, shifting it 
to the UV. 

The differences between the spectral fluxes for active compo- 
nents and for the quiet Sun integrated from 900 A to 40 000 A are 
presented in Table[JJfor the LTE (A05 and G91 abundances), for 
the NLTE with G91 abundances and NLTE with enhanced col- 
lisions and combined A05 and G91 abundances (NLTE comb.). 
We conclude that the NLTE calculations (both with and without 
ODF) significantly reduce the variations in the TSI. 

6. Conclusions 

We have presented a further development of the radiative trans- 
fer code COSI. The code accounts for the NLTE effects in sev- 
eral hundred lines, while the NLTE effects in the several million 
other lines are indirectly included via iterated opacity distribu- 
tion function. The radiative transfer is solved in spherical sym- 
metry. The main conclusions can be summarized as follows. 

- The inclusion of the molecular lines in COSI allowed us to 
reach good agreement with the SORCE observations in the 
main molecular bands (especially in the CN violet system 
and CH G band). It has also solved the previous discrepan- 
cies between the LTE calculations with the COSI code and 
ATLAS 12 calculations. We showed that their strong temper- 
ature sensitivity allows molecular lines to significantly con- 
tribute to the solar irradiance variability. 

- We introduced additional opacities into the opacity distribu- 
tion function. It allowed us to solve the well-known problem 
of overestimating the UV flux in the synthetic spectrum. We 
should emphasize, however, that the magnitude and behav- 
ior of the additional opacity strongly depend on the applied 
model. 

- We have shown that the concentration of negative hydro- 
gen is strongly affected by NLTE effects as explained in 
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Fig. 10. Upper panel: SORCE observations vs. calculations with the FALC atmosphere model and described in the text abundances. 
Lower panel: ratio between the calculations and measurements. 
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Fig. 12. The flux differences between plage and quiet Sun (upper panel), and bright network and quiet Sun (lower panel) calculated 
in NLTE and LTE. 



Sect. 14.2.11 It decreases the level of the visible and infrared 
continuum and leads to a discrepancy with the measured 
level if the calculations are done with the current models of 
the quiet Sun atmosphere. 
- We presented calculations of the total and spectral solar irra- 
diance changes due to the presence of the active regions and 
showed that NLTE effects can strongly affect both of these 
values. 
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Fig. 13. The spectra calculated with COSI (solid line) vs. SOLSPEC measuremets (dotted line). 
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